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Abstract 
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^j^We study the effect of uniaxial strain on the electronic band structure of gapped graphene. We consider two types of gapped 
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graphene, one which breaks the symmetry between the two triangular sublattices (staggered model), and another which alternates 
the bonds on the honeycomb lattice (Kekule model). In the staggered model, the effect of strains below a critical value is only 
a shift of the band gap location. In the Kekule model, as strain is increased, band gap location is initially pinned to a corner of 
the Brillouin zone while its width diminishes, and after gap closure the location of the contact point begins to shift. Analytic and 
numerical results are obtained for both the tight-binding and Dirac fermion descriptions of gapped graphene. 
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1. Introduction 

Recently, strain engineering of the electronic structure has 
been explored as an alternative method in the design of 
graphene-based electronic circuitry |1]. The approach is based 
on generating local strains to change the hopping amplitudes in 
an anisotropic way, which in turn leads to the presence of ef- 
fective gauge fields for the Dirac electrons |2]. A key finding is 
that for small and moderate uniaxial deformations the gapless 
Dirac spectrum is robust, and a gap opens only for large defor- 
mations above a particular threshold IB U, |51] . More precisely, 
the presence of anisotropy in the tight-binding hoppings on a 
honeycomb lattice makes the Dirac points approach each other 
until they merge at a critical asymmetry, at which point a band 
gap begins to open l6l|7iy,|9D. 

For practical purposes, graphene can be considered as a 
gapless semiconductor Nevertheless, gaps of various origins 
have been identified in graphene. Intrinsically, spin-orbit cou- 
)lin g is responsible for a tiny gap, on the order of 10"^ meV 
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1311 . and electron-electron interactions may ren- 
der graphene an insulator in vacuum 1IT4I1 . Extrinsically, in- 
teraction with substrates and adlayers can induce a band gap 
in graphene. Epitaxial graphene has a band gap llSll which 
has been explained as substrate-induced ITgIi . Other substrates 
studied include boron-nitride and Cu \ Ni fl^, and several 
other metals lfl9ll . Band gaps induced by the adsorption of wa- 
ter and ammonia molecules |20], and alkali metals [21] have 
been studied theoretically. Generally, there are two ways of in- 
ducing a gap in monolayer graphene. One way is the mixing of 
electronic states with different pseudospins in the same valley, 
and another way is the mixing of states that belong to different 
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valleys. The former can be achieved by sublattice symmetry 
breaking which leaves the A and B carbon atoms in different 
environments, while the latter is produced by certain transla- 
tional symmetry breakings i2lll . 

The purpose of the present work is to study the effect of uni- 
axial strain on the band structure of gapped graphene 12211 . As 
we show, the effect depends on the origin of the band gap. If it is 
induced by AB sublattice symmetry breaking, small and moder- 
ate strains do not change the width of the gap, but cause its loca- 
tions in A:-space to move in a similar way as those of the corre- 
sponding Dirac points in the gapless case. On the other hand, if 
the band gap is induced by translational symmetry breaking that 
couples the different valley states, a different behavior emerges. 
With increasing uniaxial strain, the band gap first diminishes 
and closes in its fixed location and, afterwards, the neutrality 
point is shifted in reciprocal lattice as in the gapless case. 

Our paper is organized as follows. Section |2] contains our 
analysis of the relevant tight-binding models. Section [3] con- 
tains analysis of the corresponding Dirac equations. Section |4] 
contains numerical results and their discussion, and Section |5] 
presents our conclusions. 



2. Tight-binding model 

The models we use can be described as different types of 
strain distributions ll23ll shown in Fig. [T] Sublattice symme- 
try breaking, shown in Fig. can be associated with an 
out-of-plane strain distribution. The so-called Kekule distor- 
tion, shown in Fig. [Ifb), breaks the translational symmetry in 
a way which can be described by a V3 x V3/?30° commen- 
surate lattice, and can be the result of an in-plane strain dis- 
tribution. Finally, the quinoid distortion of graphene, shown 
in Fig. [He), represents in-plane uniaxial strain parallel to the 
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Figure 1: (a) Staggered, (b) Kekule, and (c) quinoid models. 




Figure 2: (a) Primitive unit cell of graphene. The vectors 6 defined as vectors 
from atom A[ to its nearest neighbors on the B sublattice are given by Si = 
(0,-1/ V3)a, i?2 = (1/2, l/2V3)oand(?3 = (-1/2, 1/2 V3)a. (b) Kekule-type 
graphene primitive unit cell containing its six atom basis. The Kekule distortion 
is characterized by two kinds of hoppings, one set ai'ound the hexagonal ring 
(red), and the other set crossing the unit cell boundaries (blue). The addition of 
a uniaxial strain in the vertical direction results in four hoppings, ti, t2, and 
(4. (c) The relationship between the two Brillouin zones. The path TKMT is 
used in the band plots. 



nearest-neighbor bond in the vertical direction. Our tight- 
binding models based on this description are defined by on-site 
energies and the nearest-neighbor hoppings, but we neglect the 
change in bond lengths and use the perfect honeycomb lattice. 

The primitive unit cell of graphene and the nearest neighbors 
of the A and B atoms in the unit cell are shown in Fig. |2ja). 
The tight-binding models involve the vectors 6 from an A atom 
to its nearest-nearest neighbors which are defined in the cap- 
tion of Fig. 121 The primitive unit cell of the Kekule model is 
shown in Fig.|2lb). The Brillouin zone of the Kekule model and 
that of graphene, the larger hexagon, are depicted in Fig.|2jc), 
with special symmetry points defined. In the Kekule model 
the points K, K' and F belong to the reciprocal lattice and are 
therefore equivalent, which results in the coupling of the two 
inequivalent valleys that is responsible for the band gap. 

A few parameters are needed to define the tight-binding mod- 
els. In addition to the hopping parameter t -2.7 eV for 



graphene), there is the band gap of gapped graphene and the 
asymmetry in hoppings caused by uniaxial strain. We denote 
by A the half-width of the band gap, and by 6t the shift in one 
of the hoppings which is in the dkection of strain axis. The pa- 
rameter t may be taken to be positive without loss of generality, 
and it can serve as the energy scale, so that the models contain 
two adjustable parameters A/t and 6t/t. 

The gapped graphene represented by Fig. [TJa) can be de- 
scribed by the tight-binding model with staggered on-site en- 
ergies +A. Including the effect of uniaxial strain, we can write 
the Hamiltonian of this staggered model as 
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where 
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Diagonalizing the Hamiltonian defined by Eq. ([T]), we find the 
energy bands, 



e^(k) = + + \h(k)\^. 



(3) 



which show that the band gap is > 2A. The function h(k) is the 
same as in gapless graphene under strain, where zero modes 
exist for dtjt < 1 and are shifted from K and K' toward M 11241] . 
Setting = 0, and taking ^y^f = 27r/ V3 andfejO = 2n/3+pxa, 
to move along the KK' line, we find the shift from K to be 



Pjc = - cos 
a 



2k 
3a 



(4) 



We now turn to the gapped graphene based on Kekule distor- 
tion, shown in Fig.fTfb), which has zero on-site energies but two 
alternating h oppings , on one-third and two-thirds of the bonds, 
respectively 112 U 12511 . 



2^ 
t + -A, 
3 
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where, as in the staggered model, A is half of the energy gap. 
When uniaxial strain is applied to the Kekule model, we assume 
that a shift of 5t is added to the hoppings that are parallel to the 
strain axis. We define these hoppings as 



f3 = f 1 -H 6t, 



tA-t2+ 6t. 



(6) 



The Hamiltonian of the Kekule model is given by a 6 x 6 
matrix, 

„ a. { 03x3 HAB{h\ 
HsAik) 03x3 , 

where the 3x3 block Hab is 
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't2e'^-'^ 




he'^-'^' 


HABik) = 














t2e'^-'\ 
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and HBA(k) - H^g{k). The matrix elements of Eq. (O have the 
general form tij exp(ik ■ tjj) when linking A,- and Bj atoms, and 
can be read off Fig. |2lb). 
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The band structure of the Kekule model can be obtained by 
numerical calculation of the eigenvalues of Eq. (|7]). However, 
the existence of zero modes and their locations are determined 
more simply by 



det[HK(k)] = \det[HAB(k)]\ = 0. 



(9) 



Since the zero modes are expected to occur on the horizontal 
Une through F, we calculate the determinant of Hab for ky = 0, 

Dik,- A/t, 6t/t) = 2 fi fj cos + t3(tj -t\)-2 t\. (10) 



Setting D = 0, we find 

3 A; I a 



2f4f: 
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which has a solution for ky^ if the right-hand side is in the inter- 
val [-1, 1]. In particular, if 5f = A then t\ - t\ and the right- 
hand side of Eq. (fTTI) becomes unity yielding A:^ - which is 
the r point, or its equivalent K points. Therefore, the gap in the 
spectrum at K vanishes as 6t is increased from to A. Further 
increase of 6t then causes the contact points to shift away from 
the K points, until they merge at M, when the right-hand side 
of Eq. (fTTT l becomes -1 and k^^a - 27r/3. 



3. Dirac equation 

We can use the Dirac equation to describe the effect of strain 
on low-energy electrons provided that both A/f and 5tjt <K 1. 
The Dirac equation can be derived from the tight-binding model 
by setting k = near the K and K' points. For the staggered 
model under strain we obtain 
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where p - Px + ipy and we have used units such that h - I 
and vp - ta V3/2 = 1. However, we restore vp explicitly in 
some of the derived results below. Here we have followed the 
convention 

ilf = [^KA, ^KB, ^K'A,^K'bV 

for the four-dimensional spinor |26]. Our results can be ex- 
tended to the band gap due to spin-orbit interaction which can 
be described with a similar Hamiltonian as Eq. (fT2l i. except that 
the gaps have opposite signs for K and K' points |10]. The 
energy dispersions derived from Eq. (fT2l l are given by 



e^(p) = ± ^ 



+ 6t)^ + pj + A2, 



(13) 



and can be easily verified to be the limiting cases of Eq. (|3]l. 
The shifts of the Dirac points are p^ = +6t/vf, which agree 
with Eq. © for 6t/t « 1. 



For the Kekule model under strain we have 



and the energy eigenvalues are given by 
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For 6t < A there is a gap of size A — 6t p - Q. If 6t > A, zero 
modes exist at 



Py = 0, 



: ^5fi - A2. 



(16) 



For 6t - the energy dispersions, ( fT3l l and (fTsT i, give identi- 
cally gapped Dirac spectra, with a density of states (DOS) per 
unit area given by 



P(e) 



0, 



lei < A 



2|e|/7rv?., otherwise. 



(17) 



However, for 6t Q the DOS of the staggered model remains 
the same, while that of the Kekule model changes as the gap 
shrinks. For 6t - A and e <k A, Eq. ( fTsT i becomes 



2 Px , 2 



which yields a DOS given by 

2F(l/4) VA^ 



P(e) = 



;r3/2r(3/4) ' 



(18) 
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4. Numerical results and discussion 

Equation (fTTT l can be solved numerically for k^ as a function 
of A and 6t. The results for A = 0, t/2, t are shown in Fig. [3] 
For A = the same curve can be obtained from Eq. (|4|l. For 
nonzero A, the gap first diminishes at K as 6t is increased from 
to A, and with further increase of 6t the Dirac point moves 
toward M. 

In Fig.|4]we make plots of Eqs. (fT3T l and ( fISl l for a few values 
of St, using A as the scale of energy and momentum. For 6t - 0, 
shown in Fig. Ufa), both the staggered and Kekule models give 
the same gapped spectrum. However, for St - A, the valence 
and conduction bands shift laterally in the p^ direction with- 
out a change in the gap for the staggered model as in Fig. Hfb). 
In contrast, for 6t - A the gap closes in the Kekule model for 
one pair of bands while the other pair are repelled by 2A, as in 
Fig- Etc). We note that the dispersion is quadratic instead of lin- 
ear near the degeneracy point as can be expected from Eq. (fTST l. 
Further increase of 6t to 2A, shown in Fig.HJd), causes a shift 
of the neutrality points in the p^ direction, and the dispersions 
become linear where the bands cross. 

Figure |5] shows the band structures from Eq. Q for A = 
and f/10. For 6t - 0, Figs. |5|a) and (b), the Dirac point and 
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Figure 4: Dirac fermion dispersions along the p,, = line, (a) Gapless and 
gapped Dirac fermions, with gap equal to 2A. (b) Strain of r5? = A on gapless 
and staggered models, (c) Strain of 5? = A on Kekule model, (d) Strain of 
5t = 2A on Kekule model. 



Figure 5: Band structures from the tight-binding model. (a,c,e) show the band 
structures of graphene for St = 0, t/2, t, and (b,d,f) show the band structures 
for the staggered model for the same set of strains. 

the gap, respectively, are located at the K point. For St - t/2, 
Figs.|5jc) and (d), they are shifted by the same amount to some- 
where along the KM line. For 6t = t, Figs.Qe) and (f), the shift 
reaches the M point. The last cases are critical in that the Dirac 
points merge at M. We can see that the dispersions are quadratic 
near M on the KM line, but linear on the MF line, i.e., a flatten- 
ing of the Dirac cones takes place which can be seen in contour 
plots of the band structure (as shown in Ref. 6). For St > 1 a gap 
opens at the M point for graphene, and the gap of the staggered 
model becomes wider 

Figure |6] shows the band structures for the Kekule model 
obtained from the numerical evaluation of the eigenvalues of 
Eq. O. The band structures are plotted in the extended scheme 
so they can be compared with those of Fig. |5] The FPKMQF 
path is implied in Fig.|2lc), where P and Q are the points where 
the path crosses the Brillouin zone of the Kekule model. It 
must be remarked that this path does not enclose the irreducible 
wedge of the Brillouin zone for strained graphene ||5l|27tl, but 
includes the KM line where the band crossing may occur. Fig- 
ure|6ja) may be compared with Fig.|5lb) which shows the band 
structure of the staggered model for the same A. The regions 
of the gap at the K point are similar in both figures, consistent 
with the Dirac equation description, but there are extra gaps in 
the band structure of the Kekule model at P and Q. As we have 
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Figure 6: Band structures of the Kekule tight-binding model for 6t 
0, tjW, til, 1.00299f. The extended zone scheme is used. 



seen in Fig. ID the gap closes when i5f = A and this can also be 
seen in Fig. |6tb). Here the dispersion at K is linear on the PK 
line and quadratic on the KM line. A comparison of Figs. |6jc) 
and (d) with Fig.lSjc) and (e) shows that in the Kekule model, 
after the gap closes, the behavior of the neutrality points are not 
very different from that in the gapless case. 

5. Conclusions 

We considered two models of gapped graphene, denoted by 
staggered and Kekule, respectively, and studied the effect of 
strain on their band structures. We found that in the staggered 
model the width of the band gap does not change for strains less 
than a critical value, but its locations move following the mo- 
tion of the Dirac cones in the gapless case. The effect of strain 
on the band structure of the Kekule model is less trivial. With 
increasing strain, reflected in changes in the hoppings, the gap 
begins to diminish with its locations pinned to the K points. The 
gap closes when 6t - A, i.e., when the shift in hoppings due to 
strain equals the half-width of the original gap and, afterwards, 
increasing the strain makes the neutrality point to move in k 
space as in the gapless case. The effects we have discussed may 
be observed experimentally in gapped graphene on Ni substrate 
and epitaxial graphene on SiC, respectively. 
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